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We derive the large particle number limit of the Bethe equations for the ground state of the 
attractive one-dimensional Bose gas (Lieb- Liniger model) on a ring and solve it for arbitrary cou¬ 
pling. We show that the ground state of this system can be mapped to the large- IV saddle point of 
Euclidean Yang-Mills theory on a two-sphere with a U(N) gauge group, and the phase transition 
that interpolates between the homogeneous and solitonic regime is dual to the Douglas-Kazakov 
confimenent-deconfinement phase transition. 


INTRODUCTION 

The Lieb-Liniger (LL) model is an interesting labora¬ 
tory to study properties of strongly interacting quantum 
many body systems, both experimentally mm and the¬ 
oretically. The attractive version has been used exten¬ 
sively to study, e.g., quench dynamics EH3] and is known 
to undergo a phase transition at large particle number 
m- Yet at the same time, the model is integrable and 
can be solved exactly using the Bethe ansatz Kl¬ 
in practice however, a closed form expression for the 
Bethe state has only been available in the weak and 
strong coupling limits. Therefore studies of the phase 
transition have mostly resorted to mean field methods or 
numerical diagonalization of the Hamiltonian (see, e.g., 

mm)- 

In this letter, we derive the continuum limit of the 
Bethe equations for the ground state of the attractive LL 
model and solve it for arbitrary coupling. We confirm the 
second order phase transition by considering the ground 
state energy m- 

Finally we observe that the ground state can be 
mapped exactly to the large-N saddle point of U(N) 
Yang-Mills theory on a two-sphere, where the phase tran¬ 
sition manifests itself as the confinement-deconfinement 
phase transition of Douglas and Kazakov m , which is 
deeply connected to random matrix theory |18] and has 
diverse manifestations mm- 


BETHE ANSATZ FOR THE GROUND STATE 

We are interested in the Lieb-Liniger Hamiltonian ra 

# = ~X §3~ c ^2 S ( x i- x i) W 

i<N 1 ijtj 

for N identical Bosons on the interval [0, L) with peri¬ 
odic boundary conditions. We use the sign convention in 
which the coupling is attractive for c > 0. 

The eigenstates of this Hamiltonian can be computed 
exactly using the Bethe ansatz mm- They are fully 


characterized by a set of complex numbers kj, roots of 
the Bethe equations 


ki — kj — ic 
ki — kj + ic 

jjii J 



( 2 ) 


The energy and momentum eigenvalues of the Bethe 
states are E = JY kf, P = JY kj. 

In the repulsive regime, it was shown (22] that the 
Bethe states with real roots form a complete set of the 
IV-particle Hilbert space. For attractive c > 0, there ex¬ 
ist bound states: subsets of roots with identical real but 
differing imaginary part. The ground state of the sys¬ 
tem is then characterized by a bound state of zero real 
momentum, with all the kj purely imaginary [231 . 

Replacing kj —> —ikj , we then arrive at the form of 
the Bethe equations that we will deal with: 

( 3 ) 

Ay 2 Ay L 
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Their real solution characterizes the ground state of the 
system. 

In the thermodynamic limit, where N and L are taken 
to infinity with fixed density N/L, the Bethe roots form 
an exact string with exponentially small deviations: kj ss 
c(j ~ (Y + l)/2) |23j- This limit, however, is inherently a 
strong coupling limit, which can be seen by noticing that 
by using dimensionless roots in the Bethe equations, the 
coupling becomes cL. 

At weak coupling for c —> 0 (Y fixed), the distance 
between adjacents roots is much larger than c and the 
roots are distributed according to Wigner’s semi-circle 
law [24] . 

In both limits, the roots obey the following inequality 
that is a direct consequence of , but will have to be 
imposed in the continuum limit. 


\ki~kj\>c (4) 

To see why this inequality holds, imagine changing 
the coupling adiabatically. Assume that at some point 
fcj+i — kj = c for a pair of roots (if there are several 
such pairs, focus on the one with lowest index). The sum 
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FIG. 1: Continuum root distribution p = 1 /k' (black) and numerical data points at N = 400 (red fill) 


in the Bethe equation (|3j) for ki has just one diverging 
contribution at this point and can no longer be satisfied. 

LARGE-IV LIMIT 


Studying the large N limit of the exact string solution 
reveals, that the continuum root distribution may satu¬ 
rate the constraint p(k) = 1 on an interval k £ [— b,b}. 
We therefore make the following ansatz in the solitonic 
regime 


From here on, we will consider the large N limit of the 
LL model at finite effective coupling g 

N —» oo while g = cLN = const (5) 

which is the correct limit to observe the phase transition. 
We also set L = 1 for convenience. It turns out that the 
root distribution converges to a continuous function in 
this limit if we define 


p{k) 


1 k £ [—b, b] 

p{k ) k £ [—a, —b) U ( b , a] 


( 11 ) 


Inserting into © then produces an integral equation 
for y5, the solution of which is mm 


-* ,,ni (^p) (i2) 


ki = gk(i/N) (6) 


The sum on the rhs of the Bethe equation can be 
split in a near contribution from \j — i\ < eN and the 
rest, for some e > 0. A closer analysis reveals, that for 
k' > 1, the near contribution vanishes in the double limit 
lim^o liinyv->.oo (see appendix). For the rest of the sum, 
the following Taylor expansion is valid 


kj kj + n — 2g 1 
° S ki ~ N ki- kj ' 


( 7 ) 


In the continuum limit, the Bethe equation thus be¬ 
comes an integral equation 


gk = 2 V j 


1 —fcn 


P(u) 

k — u 


du 


( 8 ) 


where p{k) = 1 /k' is the density of roots and the prin¬ 
cipal value symbol is a remnant of the e excision. The 
bounds must be chosen such that J_T ax p(u) du = 1, and 
the density must satisfy the constraint Q, which in the 
continuum limit becomes 


p(k) < 1 


(9) 


The solution of the integral equation ([8]) is Wigner’s 
semi-circle law 


p{k) = 


g 2 k 2 


( 10 ) 


As long as g < 7r 2 , the constraint © is satisfied and 
(101 is the correct ground state root distribution. For 
<7 > 7r 2 , however, this distribution violates the constraint. 


and the parameters a and b are determined from the fol¬ 
lowing conditions: 

4K(x)(2E(a:) - (l-a;)K (xj) = g 
ag = 4K(a:) and x = b 2 /a‘ 2 


where E(a;) and K(x) are the elliptic functions of the first 
and second kind, and IIi(a:,y) is the elliptic function of 
the third kind [2B], defined as |27| 


III (x,y) = I 
Jo 


du 


(1 — xu 2 )^/1 — yu 2 Vl ~ 


(14) 


Note that for g —> n 2 , we have b —> 0 and p(tt 2 ) becomes 
a semi-circle and thus the root distribution changes con¬ 
tinuously at the phase transition. 

In figure [T] we show the continuum limit root distri¬ 
bution for several values of the effective coupling. The 
numerical results for N = 400, obtained directly from I©, 
are superimposed on the graphs and match very well. 


GROUND STATE ENERGY AND PHASE 
TRANSITION 

In the large -N limit, the energy per particle becomes 


k 2 = -g 2 I k 2 p(k) dk 


N 


(15) 


For the weak coupling solution (101 this expression is 


simple to evaluate. On the strong coupling side (12), the 


integral representation of III and contour integration can 
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FIG. 2: Ground state energy per particle. Numerical 
results for 400 particles are shown in brown. In green is the 
mean-field result in the strong coupling phase. The dashed 
line shows the thermodynamic limit [251 . 


be used to calculate the energy. After simplifying with 


(131, we get 


ha 2 


for g < 7r 2 

^8(a 2 + b 2 ) + g(a 2 — b 2 ) 2 ''j for g > it 2 

(16) 


By inverting (131 and expanding a, b and finally e as 
a power series in the effective coupling g around 7r 2 , we 
find 

- e = 9 + ^(9 - 7T 2 ) 2 + 0((g - 7r 2 ) 3 ), for g> tt 2 

We observe that e(g) and e'(g) are continuous at 7r 2 , 
whereas the second derivative is discontinuous, confirm¬ 
ing that it is indeed a second order phase transition. 


As a non-trivial check, we can compare (16) with the 


expression obtained in m using mean field theory, since 
in the large -N limit we expect mean field to produce the 
correct ground state energy |28| . In our conventions, the 
mean-field ground state energy is, for g > tt 2 

- Gnf = ^ (( 2 “ TO ) E ( m ) + (! - m ) K M) ( 17 ) 

with m determined from 


4E(m) K(m) = g 


(18) 


and although we have not succeeded in proving the equiv¬ 
alence of the two expressions (16) and © by algebraic 
means, we have checked their numerical equality at var¬ 
ious values of g to several dozen digits of precision. 


EQUIVALENCE TO YANG-MILLS THEORY ON 
A TWO-SPHERE 

We will briefly recap the large N limit of the U(N) 
Yang-Mills partition function in two dimensions as de¬ 
rived in HU. It will then be obvious how our ground state 
of the Lieb-Liniger model maps directly to this theory 
quantized on a sphere. 


The partition function of pure Yang-Mills theory on a 
two dimensional manifold of genus G and area A can be 
expressed as a sum over representations R of the gauge 
group |23 

Z G {A) = ^(dimi?) 2 - 2G e" AA2c,2(fl)/2iv (19) 

R 

where A is the ’t Hooft coupling. 

For the gauge group U(N ), the sum over repre¬ 
sentations can be expressed as a sum over Young 
tableaux characterized by a set of decreasing integers 
{ni,ri2, ...,7iAr}, the components of the highest weight. 

In the ’t Hooft large -N limit the representations may 
be characterized by a continuous function h: 

Nh(i/N) = -m+i-N/ 2 (20) 

and the partition function becomes 


Z g =o(A)= j Dh(x) eM~N 2 S eS {h]) 

Ses[h] = - [ [ log \h(x) - h(y)\dxdy+ 
Jo Jo 


AX 2 [\, ,2 , AA 2 
2 in 24 


( 21 ) 


Since the n,; are monotonic, it is clear that h(x) obeys 
the inequality h(x) — h(y) > x — y, so h'(x) > 1. 

The large-N saddle point approximation of (211 yields 
an integral equation for the density p(h) = dx/dh 

" P(s ) 


AX h = 2Vi 


h — s 


ds 


( 22 ) 


Clearly this integral equation with constraint is iden¬ 
tical to equation ([8]) that governs the Bethe root distri¬ 
bution in the ground state of the Lieb-Liniger model. 
The correspondence directly maps the density of Young 
tablaux boxes h to the density of Bethe roots k and the 
’t Hooft coupling A 2 to the effective LL coupling g. 

The phase transition at g = tt 2 in the Lieb-Liniger 
model appears as the confinement/deconfinement phase 
transition in the gauge theory. 

It is not yet clear, but an interesting open question, 
how physical observables of both systems can be related 
to each other. 

Note that this is not the first time a correspondence be¬ 
tween Bethe equations of an integrable one-dimensional 
system and Yang-Mills theory has been found [30] • In 
the known examples, the system was mapped into the 
moduli space of a supersymmetric gauge theory. And 
the Bethe roots played the role of the eigenvalues of the 
complex scalar in the vector multiplet. 

In our case, however, we map the Bethe roots to the 
components of highest weight of the representation of 
U(N) that dominates the saddle point of the partition 
function. 
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B 


B 


B 


1 2.00 6 1.97 11 1.98 

2 2.00 7 1.95 12 1.94 

3 1.99 8 1.94 13 1.90 

4 1.99 9 1.91 14 1.87 

5 1.98 10 1.72 


FIG. 3 &: TABLE I: Asymptotic behavior of A (N) with 
best fit parameters B for different couplings. 


NUMERICAL CHECKS 


We have performed numerical checks to validate our 
continuum results. To this end, we have solved the 
Bethe equations at various values of N and g , using 
the Levenberg-Marquadt solver provided by Mathemat- 
ica eh- In order to probe the convergence of the finite 
N root distribution to the analytic large N expression, 
we compute the mean square deviation 


N 


a ( w 'S) = jvE 


i—1 


h(g) -gk ( —,g 


(23) 


where k(x,g) is defined by numerically integrating (12 

To}. 

The results are displayed in Fig. [3] where we show 
A ( N , g) as a function of N for different values of g. We 
observe that A ( N,g ) behaves like A(g)N~ B ^. The best 
fit parameters B(g) are shown in table [TJ Based on these 
numbers, we conjecture that B = 2 at large N, and we 
notice that subleading (in N) effects seem to be stronger 
around the phase transition g = n 2 . 


CONCLUSIONS 


In this letter we have studied the attractive Lieb- 
Liniger model in the large-JV scaling limit. We have 
derived an integral equation ([8} that provides the con¬ 
tinuum form of the Bethe equations. Together with a 
bound on the root density, this has allowed us to calcu¬ 
late the limiting form of the Bethe root distribution (10 
12 }. 

The phase transition from the homogenous weak cou¬ 
pling to the bright soliton phase manifests itself in a 
change in the functional form of the root distribution. 


The ground state energy (161 coincides with the mean- 


field result - involving an identity of elliptic integrals - 
and confirms the order of the phase transition (second 
order). 

The equivalence between the large-IV saddle point of 
U(N) YM theory on a sphere with the scaling limit of 
the LL model seems like a promising avenue for future in¬ 
vestigations, especially considering that various relations 


between (supersymmetric) YM theory and integrable sys¬ 
tems have already been uncovered EDI- 

Another interesting direction is how to compute the 
lowest-lying excitations in the large -N limit. The knowl¬ 
edge of the root density for the first excitations would 
allow us to probe the time evolution of observables at 
arbitrary couplings. 
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Appendix 


For some e > 0 the near contribution to the sum ^ 
from indices \j — i\ < eN can be written in the following 
form 

E / 2 ki-s + kt+s — 2ki 

0=1 v ^ + 

We switch to continuum variables ([6]), with x = i/N 
and assume that k'{x) > 1 and k"(x ) ^ 0. Defining 
^min = m ' n k '(y ) 2 for y € [x — e,x + e], the above sum 
can be bounded by 


eN 


Ye l0 § 



2 c\k"(x)\ 
N k’l;„ - 1 


< e 


c\k"{x)\ 

k 12 — 1 


for some c > 1. Now it is obvious that this near contri¬ 
bution vanishes in the limit lim e _>.o li m JV->oo- 





